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ABSTRACT 

In this paper we report unusual glitches in two young pulsars, PSR J1825— 0935 
(B1822-09) and PSR J1835-1106. For PSR J1825-0935, a slow glitch characterised 
by a temporary decrease in the slowdown rate occurred between 2000 December 31 
to 2001 December 6. This event resulted in a permanent increase in frequency with 
fractional size Av jv ~ 31.2(2) x 10~ 9 , however little effect remained in slowdown rate. 
The glitch in PSR J1835-1106 occurred abruptly in November 2001 (MJD 52220±3) 
with Av/v ~ 14.6(4) x 10~ 9 and little or no change in the slow-down rate. A signifi- 
cant change in i) apparently occurred at the glitch with v having opposite sign for the 
pre- and post-glitch data. 
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1 INTRODUCTION 



A pulsar glitch is a phenomenon in which there is an 
abrupt increase in the rotation frequency v, often accom- 
panied by an increase in slow-down rate. Typically the frac- 
tional increase in pulsar rotation frequency is in the range 
Av/v — 10 -9 ~ 10 -6 , and the relative increment in slow- 
down rate is Av/v ~ 10~ 3 , where v and z> are pulsar rota- 
tion frequency and frequency derivative, respectively. Larger 
glitches in younger pulsars are usually followed by an expo- 
nential recovery or relaxation back toward the pre-glitch fre- 
quency, while for older pulsars and small glitches the jump 
tends to be permanent. Glitch activity, defined as the accu- 
mulated pulse frequency change due to glitches divided by 
the data span, is high in young pulsars with characteristic 
age of 10 4 - 10 5 yr, while it is low for very young and very 
old pulsars (Shemar & Lyne 1996; Wang et al. 2000). 

The trigger of the pulsar glitch is not well understood. 
In the classic starquake model, as a consequence of long-term 
slow-down in spin rate, deformation stress in the rigid crust 
builds up to resist the decreasing oblateness (Baym et al. 
1969). When the stress exceeds a critical point, the crust 
cracks suddenly, resulting in a sudden increase in spin rate. 
In the superfluid vortex unpinning and re-pinning model, 
triggering of the glitch is due to coupling of crust and the 
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superfluid interior as a consequence of a sudden unpinning 
of vortex lines and the post-glitch relaxation is due to the 
vortex gradually re-pinning to the crust lattice (Anderson 
& Itoh 1975; Alpar, Cheng & Pines 1989). Based on the 
observed typical glitches, both of the models have a sud- 
den increase in rotation frequency and slow-down rate (i.e. 
Av/v > 0) at the time of the glitch. The post-glitch re- 
laxation represents a return to equilibrium with a linear 
response of the interior superfluid, while the lack of relax- 
ation represents a nonlinear response of the superfluid (Al- 
par, Cheng & Pines 1989; Ruderman, Zhu & Chen 1998). As 
more glitches were detected, it became clear that glitch be- 
haviour varies in aspects such as glitch rate, amplitude and 
relaxation. In some cases discrete timing behaviours such as 
slow glitches were observed, i.e., the pulsar is spun- up over a 
time scale of days, weeks or even months (Shabanova 1998; 
Wang et al. 2000; Wong, Backer & Lyne 2001), accompanied 
by decreased slow-down rate (Av/v < 0). These diverse fea- 
tures suggest glitches are triggered locally in the superfluid 
interior. 



In addition to glitches, pulsars also suffer another kind 
of timing irregularity known as timing noise, which is charac- 
terised by restless, unpredictable, smaller scale fluctuations 
in spin rate (Cordes & Downs 1985) with timescales from 
days to years. The timing noise induced fluctuations of pulse 
frequency are small, with fractional changes Sv/v < 10 -9 . 
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Younger pulsars generally show more timing irregularities 
(Cordcs & Downs 1985; D'Amico ct al. 1998). 

In this paper we report two unusual glitch events in PSR 
J1825-0935 (B1822-09) and PSR J1835-1106 detected 
by Urumqi 25 m radio telescope. PSR J1825— 0935 is well 
known by its rare properties of interpulse, mode-changing, 
drifting sub-pulse and microstructure in the pulse emission 
(Fowler et al. 1981; Gil et al. 1994). Earlier observations have 
revealed unusual glitch events for PSR J1825— 0935 (Sha- 
banova 1998, Shabanova & Urama 2000), in which there was 
a gradual spin-up or equivalently a decreased slow-down rate 
over a few hundred days. In this paper we present a new but 
similar discrete timing event which occurred between 2000 
December 31 and 2001 December 6. We interpret this as a 
slow glitch as suggested by Greenstein (1979) and Cordes 
(1979). 

PSR J1835— 1106 is a young pulsar (characteristic age 
~ 10 5 yr) discovered in Parkes Southern Pulsar Survey 
(Manchester et al. 1996) which has no previously detected 
glitch. In Section 2 we introduce the observations and data 
analysis, Section 3 describes the glitches in detail, and in 
Section 4 we discuss our results. 



2 OBSERVATIONS AND DATA ANALYSIS 

Observations for PSR J1825-0935 and PSR J1835-1106 
were made regularly at 1540 MHz as part of the Urumqi Ob- 
servatory timing program which commenced in 1999. Cur- 
rently about 200 pulsars are monitored with an average in- 
terval between observations of about nine days (Wang et al. 
2001). The data for this work consist of observations span- 
ning ~1400 days from January 2000 to November 2003. The 
receiver was a dual-polarisation room-temperature system 
before July 2000 and was then upgraded to a cryogenic sys- 
tem. The dedispersion is provided by a 2 x 128 x 2.5 MHz 
filterbank. The 1-bit digitised data were sampled at 1-ms in- 
tervals and on-line folded at the topocentric period to form 
sub-integrations. The offline data reduction includes dedis- 
persing and summing the data in frequency and time to form 
mean pulse profiles. 

The mean pulse profile from each observation was 
cross-correlated with a high quality template to produce 
pulse topocentric times of arrival (TO As), which were then 

analysed using the timing program TEMPOS. The JPL 
ephemeris DE200 (Standish 1982) was used to correct TO As 
to the Solar-System barycentre. The time-corrected data 
for a given pulsar were fitted with the standard spin-down 
model in which the predicted pulse phase §(t) is expressed 
as: 

$(t) = $o + M* - to) + \i>o(t - to) 2 + \i>o{t - to) 3 , (1) 
2 6 

where <&o is the pulse phase at time to, and vo, i>o and i>o 
are the pulsar rotation frequency, frequency derivative and 
frequency second derivative at time t respectively. 

In a classic glitch, there is a sudden deviation in the ob- 
served pulse phase at the time of the glitch due to the jump 
in frequency and frequency derivative, plus an exponential 
decay in the frequency jump as a function of time. However 
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as will be described in more detail in Section 3.1, the glitch in 
PSR B1825— 0935 builds up gradually and we describe the 
glitch effect by comparing the timing solutions away from 
the spin-up event. For PSR J1835— 1106, the analysis shows 
an abrupt jump in frequency but little change in frequency 
derivative at the glitch epoch, and no exponential decay. 

Table 1 gives the pulsar names and J2000 positions in 
the first three columns, the remaining columns contain the 
pulsar dispersion measure (DM) and characteristic age. Un- 
certainties in the last quoted digit are given in parentheses. 
Given these parameters and rotation model in Equations 1 
we obtained the glitch parameters which will be described 
in detail in next section. The positions and DMs were held 
fixed in the fitting process. 



3 RESULTS 

3.1 PSR J1825-0935 (B1822 09) 

Fig. 1(a) shows the timing residuals spanning the whole data 
set after fitting for v and v. The three dashed lines, repre- 
senting MJDs 51909, 52249 and 52798 respectively, divide 
the data into four sections. Large residuals arise when fit- 
ting data across the dashed lines and a change in v between 
the first two dashed lines is evident. The rotation parame- 
ters derived from each section are given in Table 2. Fig. 1(b) 
shows the timing residuals for the whole data set with re- 
spect to the timing solution for the first section. The curva- 
ture of the residuals between the first two dashed lines (MJD 
51909—52249) reveals a gradual spin-up or, more accurately, 
a decrease in the spin-down rate for PSR J1825— 0935. 

This decrease in spin-down rate is clearly shown in 
Fig. 2, which shows the variations of rotation frequency and 
frequency derivative with respect to timing solution before 
MJD 51909. The individual values of v and v in the plot 
are derived from short fits spanning 50-100 d. The plot con- 
firms the gradual increase of pulsar frequency between MJD 
51909-52249, and another possible similar but much smaller 
event near the end of the data, as indicated by the third 
dashed line. Following Greenstein (1979), we describe these 
events as slow glitches. As shown in Fig. 2(a), for the first 
slow glitch, there was a continuous increase in frequency for 
about 300 d followed by a return to the initial state for the 
next 500 d. Associated with the spin-up process is a decreas- 
ing slow-down rate (increasing v) which lasted ~ 120 d with 
maximum Av ~ 3 x 10~ 15 s~ 2 . The slow-down rate then de- 
cayed to approximately the pre-glitch level within ~ 220 d. 
The fractional changes in frequency and frequency derivative 
before and after the slow glitch are Av fv = 31.2(2) x 10~ 9 
and Aviv — 1.9(1) x 10~ 3 respectively. These values and 
the approximate size of the second possible slow glitch are 
given in Table 3. 

3.2 PSR J1835 1106 

Fig. 3(a) shows the pre-fit residuals spanning the whole data 
set with respect to the timing solution before MJD 52220. 
The sudden change in residual slope after MJD 52220 indi- 
cates a clear glitch at this time. The residuals after fitting 
for a glitch model with jumps in frequency and frequency 
derivative at MJD 52220 as well as the pulse frequency and 
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Table 1. The parameters of PSR J1825-0935 and PSR J1835-1106. 



PSR J 


PSR B RA(J2000) 


DEC(J2000) 


DM 


Age 




(h m s) 


(d m s) 


(cm -3 pc) 


(10'V) 


1825-0935 


1822-09 18:25:30.596(6)" 


-09:35:22.8(4)° 


19.39(4) b 


2.33 


1835-1106 


18:35:18.287(2) c 


-11:06:15. l(2) c 


132.679(3) 


1.28 



a Arzoumanian ct al., 1994 b Hobbs ct al, 2003 c D'Amico et al., 1998 



Table 2. The rotation parameters for PSR J1825-0935 and PSR J1835-1106. The errors are at 2<r level. 



PSR J 


V 


V 


V 


Epoch 


Fit Span 


Residual 


No. of 






(10- 12 s- 2 ) 


(10- 24 s- 3 ) 


(MJD) 


(MJD) 


(ms) 


TOAs 


1825-0935 


1.30039948201(3) 


-0.088487(8) 




51718.0 


51549-51886 


0.73 


49 




1.3003967479(2) 


-0.08684(2) 


-9.0(7) 


52079.0 


51909-52249 


1.46 


40 




1.30039331524(3) 


-0.088657(5) 




52529.0 


52287-52769 


0.87 


44 




1.3003906000(2) 


-0.08850(7) 




52884.0 


52798-52969 


0.93 


14 


1835-1106 


6.0273144867(3) 


-0.74902(2) 


96(4) 


51909.0 


51600-52218 


1.34 


68 




6.0272702839(4) 


-0.74807(2) 


-74(4) 


52595.0 


52221-53021 


2.12 


70 



o 




MJD - 51500 



Figure 1. Timing residuals for PSR J1825— 0935 spanning the 
whole data set (a) after fitting for v and v and (b) with respect to 
the model before the slow glitch. The slow glitch occurred between 
MJD 51909-52249 as indicated by the first two dashed lines. An- 
other possible slow glitch with smaller amplitude occurred near 
the end of the data set as indicated by the third dashed line. 



its first two derivatives are shown in Fig. 3(b). Significant 
cubic terms with opposite signs are present in the pre- and 
post-glitch data, suggesting a sign change in i> at the time 
of the glitch. Fitting separately to the data sets before and 
after MJD 52220 gave the solutions in Table 2 showing this 
sign change. 

Fig. 4 presents the variations of v and i> relative to 
a fit of v and v to the pre-glitch data. This plot confirms 
the jump in v at MJD 52220 and the reversed sign of v 
before and after the glitch. There was little or no change in 
v at the time of the glitch. A coherent timing fit across the 
glitch gives the fractional jumps in frequency and frequency 
derivative at MJD 52220 listed in Table 3. 



o rj 

CO 




500 1000 1500 

MJD - 51500 

Figure 2. Variation of v and for PSR J1825-0935. (a) Fre- 
quency residuals Au relative to the pre-glitch solution, (b) Ob- 
served variations of i>. 



4 DISCUSSION 

Shabanova & Urama (2000) discussed four glitches in PSR 
J1825-0935 between 1994 September to 1999 February. 
However we interpret them as two slow glitches, with the 
first occured during MJD 49940 to 50557 (glitches 2a to 2b 
of Shabanova & Urama, 2000), and the second slow glitch 
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Figure 3. Timing residuals for PSR J1835— 1106 (a) with respect 
to the timing solution before MJD 52220 and (b) after fitting for 
a glitch at MJD 52220. 
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Figure 4. Variation of u and v across the glitch of PSR 
J1835— 1106. (a) The frequency residual Av after subtracting a 
fit of v and v to the pre-glitch data, and (b) Observed variations 
of u. 



Table 3. The glitch parameters of PSR J1825-0935 and PSR 
J1835-1106. The errors are at 2a level. 



PSR J 


Glitch epoch 


Date 


Au/u 


Au/u ] 




(MJD) 




(10-9) 


(io- 3 ) 


1825-0935 


51909-52249 


001231-011206 


31.2(2) 


1.9(1) 




52798-52969 


030608-031126 


-2.1(5) 


— 1.8(8) 


1835-1106 


52220(3) 


011107 


14.6(4) 


-1.0(2) 



beginning at 51054 but still not completed at the end of their 
data set (glitch 4 in Shabanova & Urama, 2000). Shabanova 
& Urama (2000) state that each slow glitch was preceded by 
a small glitch (glitches 1 and 3 in their paper) with fractional 
size 10~ 10 , however we did not detect such an event before 
the third slow glitch we report in this paper. In the first slow 
glitch, the continuous increase in pulsar frequency lasted for 
620 d, leaving a permanent increase in frequency with am- 
plitude 16 nHz. The second slow glitch lasted at least 120 d, 
with frequency increase > 9 nHz. The third event reported 
here lasted ~ 340 days and is similar to the previous two 
but with a much larger amplitude of 46 nHz. The possi- 
ble fourth slow glitch lasted at least 170 d with frequency 
increase > 3 nHz. These spin-up events are separated by 
1114 d, 909 d and 549 d respectively. 

It appears that PSR J1825— 0935 switches between nor- 
mal steady slow-down and intervals of decreased slow-down 
rate. Timing noise in pulsars is usually attributed to fluc- 
tuations in the interior neutron supcrfluid and its pinning 
to the neutron-star crust (Alpar, Nandkumar & Pines 1986; 
Ruderman 1991). These two phases may represent different 
states of the interior superfluid or, less likely, of the magne- 
tosphcric configuration. 

The glitch detected in PSR J1835-1106 has the normal 
abrupt increase in pulsar rotation frequency but little change 
in v at the glitch epoch, similar to small glitches detected 
in other pulsars. A more interesting feature is the apparent 
reversal of the sign of v at the time of the glitch. Again, this 
most likely originates from a change in the properties of the 
interior superfluid at the time of the glitch. 



5 SUMMARY 

In this paper we present two unusual glitches observed re- 
cently in PSR J1825-0935 and PSR J1835-1106 at Urumqi 
Astronomical Observatory. These two glitches are both dif- 
ferent to glitches in most other pulsars, demonstrating the 
great diversity of glitch behaviours. 

The main aspects for PSR J1825— 0935 glitches are as 
follows. They are slow glitches, with the pulsar frequency 
continuously increasing for several hundred days, similar to 
the glitch events reported by Shabanova & Urama (2000). 
The persistent increases in frequency result from a tem- 
porarily decreased slow-down rate lasting several hundred 
days, after which the slow-down rate returns to its stable 
value. The main event reported here results in a permanent 
frequency increase of 46 nHz, several times larger than the 
previously reported events. 

For J1835— 1106 there was a clear glitch of relative size 
~ 10~ s at MJD 52220, the first glitch observed in this pulsar. 
Unlike glitches observed in most other pulsars, there was no 
increase in slow-down rate at the time of the glitch. There 
was however a reversal in sign of v apparently associated 
with the glitch, indicating a change in the properties of the 
process(es) responsible for timing noise in this pulsar. 
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